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Abstract 

We solve the Dirac equation in one space dimension for the case of a Unear, 
Lorentz-scalar potential. This extends earlier work of Bhalerao and Ram [Am. 
J. Phys. 69 (7), 817-818 (2001)] by eliminating unnecessary constraints. The 
spectrum is shown to match smoothly to the nonrelativistic spectrum in a 
weak-coupling limit. 
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I. INTRODUCTION 



The linear potential V{x) = g\x\ is a natural choice for a confining potential in one space 
dimension. The nonrelativistic Schrodinger equation admits a nearly analytic solution for 
this potential in terms of an Airy function and the zeros of this function and its derivative. 
The Dirac equation, on the other hand, appears to be problematic for this potential. If 
V is introduced as the time component of a Lorentz two- vector, no bound-state solutions 
exist .H'H If it is introduced as a Lorentz scalar, Bhalerao and RamQ find only a very limited 
set of solutions, with no obvious correspondence to the nonrelativistic solutions. Such an 
outcome in the scalar case is unexpected because the Klein paradox is not a problem; positive 
and negative energy particles both see a confining potential. A nonrelativistic limit for the 
positive-energy solutions should reproduce the known nonrelativistic spectrum. 

This inconsistency in the scalar case can be resolved.S The solution found by Bhalerao 
and Rami turns out to be over constrained. Here we will construct a more general solution 
and show that the nonrelativistic results are recovered in an appropriate limit. 

To see that the Dirac-equation solution should match on to the nonrelativistic solution, 
consider the equation (with h = 1 = c) 

[ap + f3{m + g\x\)]i/j = Ei/j (1.1) 

in a representation wherei 

a^a = ay= (^^. /5^/3 = a^=|^J_J^. (1.2) 

As usual, let i/j = ijj = ^ ^ ^ ^^"^ decompose the matrix equation into two coupled equations 

— v' + {m + g\x\)u = Eu , 

u — {m + g\x\)v = Ev . (1.3) 

For energies E = m + e near the rest mass m, with e ^ m, and for weak coupling g -C 
m^, the second equation yields v ^ ■u'/2m. Substitution into the first equation brings 
—u" /2m + g\x\u ~ eu, which is immediately recognized as the nonrelativistic Schrodinger 
equation. 

The solution to the Schrodinger equation is obtained by noting that u" = 2'm{g\x\ — 
e)u is the differential equation for a shifted Airy function.Q This yields the normalizable 
solution u{x) = A^Ai (^(2m5f)^/'^[|a;| —e/g]j, with Af a normalization constant. Continuity 

of u and tt' at x = requires that either Ai' (^—{2mgy^'^e/g^ = [for even solutions] or 

Ai (^—{2mgy^'^e/g^ = [for odd solutions]. Let — p„ and — p'„ denote the nth zeros of Ai 

and Ai', respectively. Then the nonrelativistic eigenenergies are e'^ = p'„((7^/2m)^/^ for even 
solutions and e„ = pn{g'^ /2'my/^ for odd. The values of p„ and p'„ can be obtained from 
tables in Ref. 0. The first four of each are as follows: p„ = 2.3381, 4.0879, 5.5206, 6.7867 
and p^ = 1.0188, 3.2482, 4.8201, 6.1633. 
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We would expect the Dirac equation to yield these same results in the limit of weak 
coupling. To see that negative energy solutions do not cause any difficulties, we can use the 
methods of Coutinho, Nogami, and ToyamaB to prove the following extension of their theo- 
rem B: For a scalar potential that is everywhere nonnegative, the positive energy solutions 
have energy E > m and the negative energy solutions have E < —m. The proof depends 
on the freedom to pick as real the solutions u and v to the coupled equations (|1.3|) . Inner 
products of u and v with the terms of these equations yield 

J vudx + J u{m + V)udx = ^ J u^dx , 

J vu'dx — j v{m + V)vdx = ^ J v'^dx , (1.4) 

where g\x\ has been replaced by a generic scalar potential V and an integration by parts 
has been performed in the first term of the first equation. For positive E and V, the second 
equation implies that / vu'dx > and then the ffist equation yields E > m. Analogous 
steps for negative E yield / vu'dx < and E < —m. Thus the two parts of the spectrum 
are completely separate, and we are allowed to focus on the positive-energy solutions only. 

II. SOLUTION OF THE DIRAC EQUATION 

To solve the Dirac equation directly, we use the same representation as Bhalerao and 
RamJi that is 



« = ^s/ = L- n ' = = T n • (2-1) 





For = ( ^ 1 they obtain the coupled equations 



V 



u + {m + g\x\)u = Ev , 
-v' + {m + g\x\)v = Eu . (2.2) 

These equations decouple in terms of their variable ^ = ^/g{m/g + \x\), such that for x > 0, 

/ d^ ^ 



'de 



+ eju = {E'/g + l)u, 
+ Av = {E'/g-l)v, (2.3) 



and for a; < 0, 



+ey={E'/g-l)u, 

+ e]v = {E^/g + l)v. (2.4) 



de 
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Obviously these are harmonic-oscillator-type equations. The normalizable solutions can be 
constructed from the Hermite functions! of order v and z/ + 1, where E'^ = 2(z/ + l)(yf, aJl 

_ rCe-«V2i7^+i(0, a;>0 
_ f Cj=e-^'I^H,{i) , X > 

However, because ^ is always positive, z/ is not restricted to being an integer.i 

Continuity at a; = requires that CHy^i{a) = C EHy{a) / ^ and C Hy^i{a) = 
CEH^{a)/ where a = mj ^ is defined as in Ref. f^. These conditions yield C = ±C 
and 

= ±^H,{a) , (2.6) 

the latter being the eigenvalue condition. This condition has a rich set of solutions when 
free of the restriction to integer u; there are infinitely many solutions for any positive value 
of a. 

The sign that appears in the eigenvalue condition ( |2.6| ) corresponds to the parity of the 
solution. The parity operatoilli is reflection in x combined with multiplication by the Dirac 
matrix (3. Because ^ is independent of the sign of x, we find that 

/ u{-x) \ _ ( v{-x) \ _,( u{x) \ , . 

^[v{-x) J ~ [u{-x) J ~^[v{x) J ■ ^ 

The presence of such a symmetry is, of course, necessary for the match to the nonrelativistic 
solution. 



III. NONRELATIVISTIC LIMIT 



To recover the nonrelativistic solution, we must take an appropriate limit. We write 
E = m(l + e) and consider small e as well as small g. The latter corresponds to large a and 
large u. In the limit of large u, we find from Ref. B that has the asymptotic form 



1/4 



Ai(t.) 



(3.1) 



where z^, = ^l\/1v + 1 and, for 



cos z, 



2/3 



(3.2) 



This immediately looks promising because the desired Airy function is present. We next 
expand in powers of and e, with yj^ = e — ^\x\/a ± 1/2q;^ and use of cos"^(l — y) 

^/2y + to obtain 
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1 + 

^ ^ ^' -i-y-, (3.3) 



1 + 2e + £2 - i/«2 

^ - 1 , Wi ^ - ( ^ 1 Z/+ . (3.4) 

Note that ^Jg\x\la is of order e at the classical turning point, which sets a natural scale for 
and that the nonrelativistic correspondence implies that e is of order {g'^/m^Y/^ ~ a~^^^. 
Thus the terms can be dropped relative to a^^. 
At lowest order, these expansions imply 

e-«'/2i/.(0, e~^'"H,^,{0 ~ Ai ({2m' / g'Yl'[\x\/ a - e]) . (3.5) 



To compare with the original nonrelativistic reduction we must connect the two representa- 
tions of the Dirac matrices. They are related by a unitary transformation 

such that 

v'=(:)=f/v=c(:). (3.7) 

Therefore, we have u = ^{u + v) and v = -^{v — u), with u and v given by (P75|). Thus for 
large a, u does indeed reduce to an Airy function with the correct argument, given e = e/m 
and a = m/ 

At the next order we obtain, with the aid of Stirling's formula! for V , 

H^iO ~ 2-VVV2r (^^) (^Jp^ (1 + ^_/4)Ai {-{2m'/g^f/^_) (3.8) 

and 

^f.,,(0~2-3/VV2^r(^) (^y'"(l + W4)Ai(-(2mV/)^/V) . (3.9) 
The combination that appears in the eigenvalue condition 



= H,+i{a) T V2z/ + 2H,{a) 

then reduces to 

^ v^2-VVV2r f ^) { ^) (3.10) 



X 



2 n , 

1 + y_/4)Ai (-(2mV/)'/'y+) T (1 + y+/4)Ai (-(2mV/)i/='y^) 



x=0 



5 



First-order Taylor expansions of the Airy functions about —{^ra^jg^^^l'^e yield 



Ai 



(3.11) 



Ai (-(2m V^?')'/^ + (2a)-2/3Ai' [-{2m^ j g")^!^ 



where terms of order higher than a have been dropped. For even parity (the upper sign) 

we have Ai' [—(^m^jg^Y^'^^ — ^^'^ odd parity, Ai [—{2m^/g'^y^^e^ ^ 0, which are 
the nonrelativistic eigenvalue conditions. Therefore the eigenvalues will match in the limit 
of large a and small e. 



IV. RESULTS AND CONCLUSIONS 



Comparison of the relativistic and nonrelativistic eigenvalues is made in Fig. where 
each is plotted as a function of l/a for the four lowest levels of each parity. The nonrelativis- 
tic values were already obtained above as explicit functions of a which are simply plotted 
as lines in the figure. The relativistic values were computed numerically, with use of Math- 
em atic A to solve the eigenvalue condition ( p.6|) at selected values of a. The dimensionless 
e is related to z/ by e = a/2z/ + 2/a — 1. For 1/a near 0, i.e. large a, the relativistic and 
nonrelativistic results are indistinguishable. As a decreases, they separate smoothly. 




From these results we see that the Dirac equation with a scalar linear potential is a well- 
defined problem in one dimension, with a rich set of solutions and a smooth nonrelativistic 
limit. As an exercise, one could extend this work to include calculation of the relativistic 
wave functions and make direct comparisons with the nonrelativistic Airy functions. They 
will match in the large-a limit. A second interesting exercise is a comparison of the ultrarel- 
ativistic, strong-coupling limit of the spectrum to the a~^^^ behavior of the nonrelativistic 
spectrum. The plots in Fig. |l| appear to imply a behavior for small a. 
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